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Abstract 

Evolution of the elliptic flow of hadrons in heavy-ion collisions at RHIC energies is 
studied within the microscopic quark-gluon string model. The elliptic flow is shown 
to have a multi-component structure caused by (i) rescattering and (ii) absorption 
processes in spatially asymmetric medium. Together with different freeze-out dy- 
namics of mesons and baryons, these processes lead to the following trend in the flow 
formation: the later the mesons are frozen, the weaker their elliptic flow, whereas 
baryon fraction develops stronger elliptic flow during the late stages of the fireball 
evolution. Comparison with the PHOBOS data demonstrates the model ability to 
reproduce the v^ 1 ^) signal in different centrality bins. 
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1 Introduction 



Elliptic flow is defined as the second harmonic coefficient V2 of an azimuthal 
Fourier expansion of the particle invariant distribution [1] 

d 3 N 1 d 2 N 

£ W = -^ [1 + 2 " lCOsW, + 2l ' 2COs( ^ , + --- 1 ' (1) 

where <fi is the azimuthal angle between the transverse momentum of the 
particle and the reaction plane, and pt and y is the transverse momentum and 
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the rapidity, respectively. The first harmonic coefficient v\ is called directed 
flow. It can be presented as 



V\ = (cos0) 




(2) 



while the t>2, which measures the eccentricity of the particle distribution in 
the momentum space, is 



(Note, that in the coordinate system applied the z-axis is directed along the 
beam, and the impact parameter axis is la beled as x-axis. Therefore, transverse 
momentum of a particle is simply p t = Jj>\ + pi ). 



The overlapping area of two nuclei colliding with non-zero impact parameter b 
has a characteristic almond shape [2] in the transverse plane. The fireball tries 
to restore spherical shape, provided the thermalization sets in rapidly and the 
hydrodynamic description is appropriate [3,4,5,6]. When it becomes spherical, 
apparently, the elliptic flow stops to develop. Therefore, v 2 can carry important 
information about the earlier phase of ultrarelativistic heavy-ion collisions, 
equation of state (EOS) of hot and dense hadronic (or rather partonic) matter, 
and is expected to be a useful tool to probe the formation and hadronization 
of the quark-gluon plasma (QGP) [7]. 

Elliptic flow of charged particles was among the first signals measured at Rela- 
tivistic Heavy Ion Collider (RHIC) in Brookhaven [8], and now there is plenty 
of data concerning the centrality, the (pseudo)rapidity, and, especially, the 
transverse momentum dependence of the v 2 in gold-gold collisions at y/s = 130 
AGeV [9] and at y/s = 200 AGeV [10]. Microscopic models based on string phe- 
nomenology and transport theory are able to reproduce, at least qualitatively, 
many features of the v 2 at ultrarelativistic energies [11,12,13,14,15,16,17,18], 
however, the quantitative agreement with the data is often not so good. Par- 
ticularly, magnitude of the distributions ^2(77) or 1*2 (pt > 1.5GeV/c) appears 
to be too high. Does it mean that the effective EOS of hot and dense partonic- 
hadronic matter in microscopic models is too soft? What should be done to 
match the data? Before starting our study it is worth noting that the ellip- 
tic flow increases due to two main processes, namely, (i) elastic or inelastic 
rescattering of particles, and (ii) non-homogeneous absorption of particles in 
spatially asymmetric dense medium. The first reaction is responsible for the 
formation and development of the hydrodynamic flow, while the latter one 
should be considered as a non-hydro contribution. 




(3) 
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Then, the microscopic calculations [19,20,21,22] show the absence of sharp 
freeze-out of particles in relativistic heavy-ion collisions. In contrast to as- 
sumptions of hydrodynamic model [23], the expanding fireball in microscopic 
models can be rather treated as a core consisting of still interacting hadrons, 
and a halo, which contains particles already decoupled from the system. The 
order of the freeze-out of different species seems to be the same for energies 
ranging from AGS to RHIC: 1 - pions, 2 - kaons, 3 - lambdas, 4 - nucleons. 
What are the consequences of the continuous freeze-out for the v 2 of these 
particles? What is the role of the rescattering and absorption in the flow de- 
velopment? When (and where) the elliptic flow is formed? 



2 Model 



To answer these questions we employ the quark-gluon string model (QGSM) 
[25,26], which is a microscopic model based on the Gribov-Regge theory [24] of 
hadronic and nuclear interactions at high energies. The main advantage of the 
GRT is the fulfillment of unitarity conditions in s- and t-channel for multipar- 
ticle processes. On one hand, the Regge-pole exchange in t-channel determines 
two-particle amplitude at s — > oo. On the other hand, Regge-pole exchanges 
fully determine multiparticle processes in s-channel, since each Regge-pole 
exchange in t-channel corresponds to a jet of hadrons with small transverse 
momenta. Therefore, unitarity is fulfilled. In the QGSM the subprocesses with 
quark annihilation and quark exchange correspond to the so-called Reggeon 
exchange in two-particle amplitudes in the GRT, while the subprocesses with 
colour exchange are represented by the one and more Pomeron exchanges in 
elastic amplitudes. This scheme is closely related to that given by the partonic 
model, where most of calculations imply somehow calculation of contributions 
of diagrams arising in the GRT. More information concerning the link be- 
tween the GRT and quantum chromodynamics (QCD) can be found, e.g., in 
[24,27,28,29] and references therein. 

Note also, that number of scattering partons (partonic density) within hadrons 
or nuclei, colliding at ultrarelativistic energies, increases with the rise of center- 
of-mass energy y/s as s a ^- ap with op(0) > 1 being the intercept of a Pomer- 
anchuk pole. This leads to the transition to Froissart regime corresponding to 
multiparton collision of black disks filled out by slow partons, i.e., the state 
which is a precursor of the Color Glass Condensate (see, e.g., [30] and refer- 
ences therein). The QGSM incorporates also the string fragmentation, reso- 
nance formation, and hadronic rescattering. The latter implies that the decay 
products of a string, - stable hadrons and their resonances, - can further in- 
teract with other hadrons. Due to uncertainty principle, secondary hadrons 
are allowed to interact again only after certain formation time, but hadrons 
containing the valence quarks of the colliding hadrons/nuclei can interact im- 
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mediately with the reduced cross section a q N taken from the additive quark 
model. Angular and momentum distributions of secondaries are adjusted to 
available experimental data. In addition, the one-pion exchange model, de- 
tailed balance constructions, and isospin symmetry arguments are used in 
case of lacking the experimental information. The positions and momenta of 
nucleons inside the colliding nuclei are Monte Carlo generated according to 
the Woods-Saxon density distribution and the Fermi momentum distribution, 
respectively. The Pauli blocking is taken into account by excluding the scat- 
tering into occupied final states. 

The transverse motion of hadrons in the QGSM arises from different sources: 
(i) primordial transverse momentum of the constituent quarks, (ii) transverse 
momentum of (di)quark-anti(di) quark pairs acquired at string breakup, (iii) 
the transverse Fermi motion of nucleons in colliding nuclei, and (iv) rescatter- 
ing of secondaries. Parameters of the first two sources are fixed by comparison 
with hadronic data. The Fermi motion changes the effective transverse distri- 
bution of strings formed by the valence quarks and diquarks of the target and 
projectile nucleons. Thus, the original strings are not completely parallel to 
the beam axis. Further details of the QGSM can be found elsewhere [25,26]. 



3 Development of elliptic flow 

To investigate the development of the elliptic flow ca. 20- 10 3 gold-gold colli- 
sions with the impact parameter b = 8fm were generated at ^/s = 130 AGeV. 
According to previous studies [8,11] the elliptic flow of charged particles is 
close to its maximum at this impact parameter, and the multiplicity of secon- 
daries is still quite high. The time evolutions of the v 2 of pions and nucleons 
as functions of rapidity 



are displayed in Fig. 1(a). Here the snapshots of the v 2 profile are taken at 2 
fm/c, 4 fm/c, ... 10 fm/c, and at the time of thermal freeze-out. This means 
that the hadronic composition was frozen at certain time t — ti, when all 
interactions were switched off and particles were propagated freely. To avoid 
ambiguities, resonances were allowed to decay according to their branching 
ratios. Surprisingly, at t — 2 fm/c elliptic flow of pions is weak. The flow con- 
tinuously increases and reaches its maximum value v\ (y — 0) ~ 5% already at 
t — 6 fm/c. From this time the elliptic flow does not increase anymore. Instead, 
it becomes broader and develops a two-hump structure with a relatively weak 
dip at midrapidity. The flow seems to continue development till the late stages 




d 2 N 
dydp t 



dp t 



(4) 
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Fig. 1. Gold-gold collisions with b = 8fm at -^/s = 130 AGeV: (a) Evolution of 
elliptic flow of pions (snapshots are made at 2, 4, ... 80 fm/c), and (b) contributions 
to the resulting (y) coming from pions frozen at 2, 4, ... 20 fm/c, respectively. 

of the system evolution. Bearing in mind the absence of sharp freeze-out of 
particles in model calculations, it is important to study the contribution of 
the survived particles to the resulting elliptic flow. The corresponding rapid- 
ity distributions presented in Fig. 1(b) reveal the peculiar feature: the vi of 
pions, which are frozen already at t — 2 fm/c, is the strongest among the 
fractions of the flow carried by pions decoupled from the fireball later on. The 
later the pions are frozen, the weaker their flow. The two-hump structure of 
the signal develops here as well, but the widths of all rapidity distributions 
are the same. For pions frozen after t = 6 fm/c the dip at midrapidity is seen 
quite distinctly. From here one can conclude that the strong elliptic anisotropy 
of pions, which left the system early, is caused by the absorption of the pion 
component in the squeeze-out direction. 

For nucleons the evolution picture of the V2(y), shown in Fig. 2(a), is similar 
to that for pions. The flow is quite weak at t = 2 fm/c, then it increases and 
gets a full strength at midrapidity between 8 fm/c and 10 fm/c, i.e. later than 
the elliptic flow of pions. Similarly to v^d/), it develops a two- hump structure, 
but the humps tend to dissolve at late stages of system evolution. In contrast 
to this behavior, the freeze-out decomposition picture of v% (y), presented in 
Fig. 2(b), does not show monotonic tendency within first 8 fm/c of the reaction: 
The flow of nucleons frozen at 2 fm/c is identical to that of nucleons frozen 
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: N evolution at 




Fig. 2. The same as Fig.l but for the elliptic flow of nucleons. 

at 8fm/c, whereas nucleons decoupled from the system between 2fm/c and 
8fm/c almost do not contribute to the resulting elliptic flow. Nucleons, which 
are decoupled after 8 fm/ c, have significant anisotropy in the momentum space, 
and the time development of flow in the nucleon sector is opposite to that in 
the pion one. Namely, the later the nucleons are frozen, the stronger their 
elliptic flow. 

To clarify the role of particle freeze-out for the formation and evolution of their 
elliptic flow, the dN/dt distributions of n c h, n, N, and A, which are decoupled 
from the system after the last elastic or inelastic collision, are depicted in 
Fig. 3(a). Here two features should be mentioned. Firstly, a substantial part 
of hadrons leave the fireball immediately after their production within the first 
two fm/c, in stark contrast with heavy-ion reactions at AGS [19] and at SPS 
[20] energies. Secondly, compared to lower energies, the particle distributions 
at RHIC are peaking at earlier times. For instance, the pion distribution has 
maximum at t^ a2 . = 5-6 fm/c, while the distributions of baryons are broader 
due to large number of rescatterings shifting their maxima to later times, 
tg ax = 10-12 fm/c (cf. « 10 fm/c and tg ax « 18 fm/c at both AGS and 
SPS energies). 

Elliptic flow carried by these hadronic species is presented in Fig. 3(b). The 
baryonic and mesonic components are completely different: pions emitted 
from the surface of the expanding fireball within the first few fm/c carry 
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Fig. 3. (a) dN/dt distribution of n^, tt, N, A over the time of their last interac- 
tion, and (b) elliptic flow of these particles for Au+Au collisions with b = 8fm at 
= 130 AGeV. 

the strongest flow, while later on the flow of pions is significantly reduced. In 
contrast to pions, the baryon fraction acquires stronger elliptic flow during 
the subsequent rescatterings, thus developing the hydro-like flow. The satura- 
tion of the flow at the late stages can be explained by the lack of rescattering 
as the expanding system becomes more dilute, and by the restoration of the 
symmetry of particle momentum distributions in the transverse plane. 

The flow has a multicomponent structure. In addition to hydrodynamic flow 
there are other components of the t>2 caused mainly by the particle splash 
from the surface area during the initial phase of nuclear collision, and by the 
non-uniform absorption of hadrons in spatially asymmetric dense matter. For 
instance, the anisotropic absorption of jets in x- and y-directions will con- 
tribute to the development of the elliptic flow of particles [31] with high trans- 
verse momenta. Hence, one may expect that the resulting elliptic flows of e.g. 
baryons and mesons obtained after the particle freeze-out are different. Figure 
4 shows the rapidity dependence of the v 2 of most abundant particle species in 
Au+Au collisions at RHIC energies, namely pions, kaons, nucleons, and lamb- 
das. Surprisingly, these distributions are almost indistinguishable: the widths 
and magnitudes of the signals are the same. To check the validity of the model 
calculations one has to perform more complex analysis of experimental data 
e.g. by studying the v^(y) distributions in different pr-intervals. Other model 
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Fig. 4. Rapidity distribution of the resulting elliptic flow of nucleons (full line), 
pions (dashed line), kaons (dotted line), and lambdas (dash-dotted line) in Au+Au 
collisions with b = 8 fm at y/s = 130 AGeV. 



predictions for the rapidity distribution of the elliptic flow of hadrons at both 
RHIC energies, 130 AGeV and 200 AGeV, can be found in [11,32,33,34]. Re- 
cently, PHOBOS collaboration presented the data on the ^(f?) distributions 
of charged particles, produced in Au+Au collisions at 200 AGeV, in different 
centrality bins [35]. Model calculations are plotted onto the experimental re- 
sults in Fig. 5. Here the elliptic flow of charged particles for the collisions with 
centralities 0-15% (0 fm < b < 2.3 fm), 15-25% (2.3 fm <b< 6.5 fm), and 
25-50% (6.5 fm < b < 9.2 fm) is presented together with the resulting flow in 
minimum bias events. We see that the model reproduces the measured signal 
pretty well, including the nearly flat distribution at \r]\ < 2 and quick fall of 
the elliptic flow at \r]\ > 2. The only discrepancy appears to arise in semipe- 
ripheral collisions at midrapidity range, where the data are approximately 15% 
above the model results. This problem can be fixed by (i) the fine tuning of 
model parameters, and (ii) proper analysis of the experimental procedure of 
data handling (for the analysis and comparison of several methods applied to 
restore the anisotropic flow see [36]). 
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Fig. 5. V2(rj) distribution of charged particles in Au+Au collisions at y/s = 200 
AGeV for (a) a/a geo = 0-15%, (b) a/a geo = 15-25%, (c) a/a geo = 25-50%, and 
(d) minimum bias events. Full symbols represent data [35], open symbols denote 
model calculations. 

4 Conclusions 

In summary, the features of the formation and development of elliptic flow in 
the microscopic quark-gluon string model can be stated as follows. First of 
all, there is no one-to-one correspondence between the apparent elliptic flow 
and the contribution to the final flow coming from the "survived" fraction 
of particles. For instance, apparent elliptic flow of pions at t — 2 fm/c is 
weak, but pions which are already decoupled from the system at this moment 
have the strongest elliptic anisotropy caused by the absorption of the pion 
component in the squeeze-out direction. Elliptic flow of hadrons is formed 
not only during the first few fm/c, but also during the whole evolution of the 
system because of continuous freeze-out of particles. Secondly, time evolutions 
of the mesonic flow and baryonic flow are quite different. Baryons, except of 
the fraction emitted within first two fm/c after the collision, are getting the 
stronger v 2 in the course of secondary interactions. Pions in average experience 
much less elastic collisions per particle because of their instant absorption 
and production in dense medium. As the substance becomes more dilute, 
their elliptic flow decreases. Freeze-out dynamics for baryons and mesons is 
also different, therefore, development of particle collective flow should not be 
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studied independently of the freeze-out picture. The general trend in particle 
flow formation in microscopic models at ultrarelativistic energies is that the 
earlier mesons are frozen, the weaker their elliptic flow. In contrast, baryons 
frozen at the end of the system evolution have stronger t> 2 . 
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